function fit_ = FitLoadingRate(p, xdata,ydata )

    %p=[R/Gamma,Gamma,dc];
    ft_ = fittype('p1*(1-exp(-p2*x))+p3',...
    'dependent',{'y'},'independent',{'x'},...
    'coefficients',{'p1', 'p2', 'p3'});

    cf_ = fit(xdata',ydata',ft_,'Startpoint',p);
        %FittedCurve = loading_rate_function(p,xdata);      
        %ErrorVector = FittedCurve' - ydata;
    fit_ =cf_;% sum(ErrorVector .^ 2);
    
    
        
end